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Abstract 

Influence of an external magnetic field on the nematic-isotropic (N — I) phase tran- 
sition in a dispersion model of nematic liquid crystals, where the molecules are either 
perfectly uniaxial or biaxial (board-like), has been studied by Monte Carlo simulation. 
Using multiple histogram reweighting technique and finite size scaling analysis the or- 
der of the phase transition, the transition temperature at the thermodynamic limit 
and the stability limit of the isotropic phase below the transition temperature for dif- 
ferent magnetic field strengths have been determined. The magnetic field dependence 
of the shift in TV — / transition temperature is observed to be more rapid than that 
predicted by the standard Landau-de Gennes and Maier-Saupe mean field theories. 
We have shown that for a given field strength the shift in the transition temperature 
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is higher for the biaxial molecules in comparison with the uniaxial case. The study 
shows that the N — I transition for the biaxial molecules is weaker than the well known 
weak first order N — I transition for the uniaxial molecules and the presence of the 
external magnetic field (up to a certain critical value) makes the transition much more 
weaker for both the systems. The estimate of the critical magnetic field (~ 110 T) for 
the common nematics is found to be smaller than the earlier estimates. 

1 INTRODUCTION 

The study of the behaviour of liquid crystals in presence of an external field (either electric 
or magnetic) has long been an active area of research because of its fundamental and techno- 
logical importance. There have been several theoretical [TJ [H HJ El (6J [7] and experimental 
[HI [9J [10l [TT1 [121 E3] reports dealing with different phase transition phenomena in presence 
of external electric and magnetic fields. 

In absence of any external field nematic liquid crystals composed of anisotropic molecules, 
usually modeled as rod-like in shape, possess long-range orientational order due to suffi- 
ciently strong intermolecular forces and in the uniaxial nematic phase (N) the molecules 
become aligned, on the average, along a single macroscopic direction called the director (n). 
The degree of this orientational order is characterized by the second rank tensor nematic 
order parameter Q [14J. With increasing temperature the order parameter of a thermotropic 
liquid crystal in the nematic phase decreases and jumps to zero as the liquid crystal under- 
goes a weakly first-order nematic-isotropic (N — I) transition at a temperature T/vj. The 
orientational order no longer exists in the isotropic phase (J). If an external field is applied, 
the director in the N phase, tends to align parallel to the direction of the field provided 
the liquid crystal material has a positive dielectric or diamagnetic anisotropy and becomes 
pinned at relatively weak fields although the orientational order increases only by a small 
amount. The transition temperature T^i also gets shifted towards a higher value. In the 
presence of an external field a weak orientational order is induced even in the isotropic phase 
resulting in an anisotropic phase known as the paranematic phase (pN). From theoretical 
studies [21 H] it has been observed that when the external field is weak the iV — I (or pN) 
phase transition occurs with a finite jump in the nematic order parameter and the transition 
is first order in nature. As the field strength is increased the difference between the N phase 
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and the pN phase decreases and finally it vanishes at a classical critical point defined by a 
critical temperature (Tc) and a critical field (-Be)- This critical end point is analogous to the 
classical critical points observed in liquid- vapour and ferromagnetic systems [15J. Above the 
critical point the I (or pN) and the N phases become indistinguishable, the symmetries of 
the phases being the same. The effect of an external magnetic field on the N — I transition 
is analogous to the effect observed by the application of pressure at the liquid-gas transition, 
the paranematic phase being the analog of the liquid-gas coexistence phase. 

There are several reports [HI QUI [HI EES] of experimentally observed electric field induced 
first-order isotropic-nematic phase transition since its earliest theoretical prediction based 
on the mean field theories of Maier-Saupe and Landau de-Gennes |14j . On the other hand 
similar studies in the presence of a magnetic field are relatively fewer. The influence of a 
magnetic field on the N — I transition was first studied by Rosenblatt [12] . His experimental 
observation demonstrated that the presence of a magnetic field of strength 14.8 T shifts the 
N — I transition temperature only by a few milliKelvin. The reason behind this small change 
is the very low value of the anisotropy of diamagnetic susceptibility for the traditional liquid 
crystal materials [H] . In order to substantially alter the N — I transition temperature for 
the conventional calamitic liquid crystals the necessary magnetic field strength should be 
larger than 100 T and the estimated critical value of the magnetic field is much higher than 

this [num. 

The experimental observation of a magnetic field induced first-order isotropic-nematic 
transition in a thermotropic liquid crystal has been reported only recently by Ostapenko et 
al [T3] who also observed a much higher shift in the N — I transition temperature (~ 0.8°C 
at B = 23 T). Observation of such pronounced effects became possible by the use of a 
new class (non-traditional) of liquid crystal molecules having a bent core and by using a 
high- field resistive magnet. However, Ostapenko et. al. could not reach the critical point 
of the system in their investigation where a magnetic field up to 31 T was used. 

The bent-core compounds have exhibited different fascinating phenomena in liquid crys- 
tal science, for example, the formation of the long searched biaxial nematic phase in ther- 
motropic liquid crystals JTSJ [Tj5] and the formation of chiral phases resulting from achiral 
molecules [20]. Another experimental study [21] has shown that the N — I transition in the 
bent-core compounds is more weakly first order as the value of T^i — T_ rj 0.4° C which 
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is significantly below the values ~ 1° C in typical calamities, T_ being the super cooling 
limit of the isotropic phase. The exceptional behaviour of this new class of liquid crystals 
has been explained [21] by considering the isotropic phase to be composed of microscopic 
complexes or "clusters" of bent core molecules. Another important aspect of the bent-core 
mesogenic molecules is their inherent biaxiality. A mean-field study [7] was reported more 
than two decades ago in which the effect of an external magnetic field on biaxial nemato- 
genic molecules was examined. It is evident from this work that as the degree of biaxiality 
is increased the critical field strength drops rapidly and for a given field strength the shift in 
TV/ is higher in uniaxial nematics composed of biaxial molecules. So far as our knowledge 
is concerned no other theoretical or computer simulation study on systems composed of 
biaxial molecules coupled with external magnetic field has been reported. 

In this paper we have tried to investigate using Monte Carlo (MC) simulation the effect 
of an external magnetic field on the N — I (or iV — pN) transition. Using a lattice model 
of nematics where the particles interact via a dispersion potential we have performed the 
multiple histogram reweighting [22] and finite-size scaling analysis [2Uj [21] of our data. We 
have studied the magnetic field dependence of the transition temperature T^i for systems 
composed of either uniaxial or biaxial molecules and have compared our results with that 
predicted by the mean field theory. We have also confirmed the first orderedness of the 
iV — pN phase transition and find that the transition gets weaker with the increase in the 
magnetic field and ultimately reaches a critical point. Using the technique similar to that 
applied by Saito et. al. in a recent work [22] for the study of quark mass dependence of 
the finite temperature QCD phase transition, we have estimated the value of the critical 
magnetic field for the system composed of uniaxial molecules. The use of finite size scaling 
along with the multiple histogram reweighting has also enabled us to obtain the magnetic 
field dependence of the supercooling limit of the isotropic phase (or pN phase ) in the 
thermodynamic limit. Expectedly we find that the difference T^j(B) — T_(B) decreases 
with the increase in the magnetic field and vanishes at a certain point (Be, Tc) which we 
have identified as the critical point. Our work is intended to supplement the predictions of 
the mean field theory in this area where no other theoretical work has been reported. The 
MC simulations which have so far been performed to study the effects of external fields on 
the N — I transitions are the works of Luckhurst and coworkers |26| and that of Berardi et. 
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al. [27]. However these simulations do not address the features of the phase transitions we 
have described. 

In the next section we describe the model employed in our simulation. This is followed 
by the computational details, results and conclusions. 

2 THE DISPERSION POTENTIAL AND THE EF- 
FECT OF THE MAGNETIC FIELD 

We consider a system of biaxial prolate molecules possessing D 2 h symmetry (board-like), 
whose centres of mass are associated with a simple-cubic lattice and subjected to an external 
magnetic field. The total energy of the system is the sum of two terms: (i) a dispersion 
potential term which takes into account the interaction between all nearest neighbour pairs 
of molecules and (ii) a field term which represents the interaction of each molecule with 
the external field. The dispersion term in the Hamiltonian contains a factor A which is 
a measure of the molecular biaxiality. The case of uniaxial molecules (D^ symmetry) is 
obtained by simply setting this parameter equal to zero. 

The orientationally anisotropic dispersion pair interaction obtained from London disper- 
sion model j28j 12H] explicitly depends on both the mutual orientation of the two interacting 
molecules (say the i th and j th molecules), and on their orientations with respect to the inter- 
molecular unit vector (r^). By isotropically averaging over the intermolecular unit vector, 
Tij the dispersion potential between two identical neighbouring molecules becomes [30| 131] 

= -%Ko(%) + 2A[i&(fy) + + 4A 2 i? 2 2 2 (^,)}. (1) 

Here fly = {<pij, %, ipij} denotes the triplet of Euler angles defining the relative orientation 
of i th and j th molecules; we have used the convention used by Rose [32] in defining the Euler 
angles, is the strength parameter which is assumed to be a positive constant (e) when 
the particles i and j are nearest neighbours and zero otherwise. R^n are combinations of 
symmetry-adapted (D 2 h) Wigner functions 

;4 = ~cos 2 0-i ( 2 ) 
Kg 

Rl 2 = ^ sin 2 ^ cos 2^ (3) 
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R 2 2Q = 1— sin 2 9 cos 20 (4) 

I i 

i?22 = j(l + cos2 0) cos 20 cos 2-0 — - cos 6 1 sin 20 sin 2-0. (5) 

The parameter A is a measure of the molecular biaxiality and it depends on the molecular 
properties. For the dispersion interactions, A can be expressed in terms of the eigenvalues 
(pi, pi, pa) of the polarizability tensor p of the biaxial molecule 

A = Jl P2 ~ P1 . (6) 
V 2 2p 3 - P2 ~ Pi ^ ' 

The condition for the maximum biaxiality is p 3 — p 2 = P2 _ Pi > 0, A = Ac = l/v6 

and this self-dual geometry corresponds to the Landau point in the phase diagram where 

a direct biaxial nematic to isotropic phase transition occurs. A < Ac corresponds to the 

case of prolate molecules whereas A > Ac corresponds to oblate molecules. This dispersion 

model can successfully reproduce both the uniaxial and the biaxial orientational orders and 

various order- disorder transitions as a function of temperature and molecular biaxiality. 

The phase diagram for this dispersion model has been well studied by mean field theory 

as well as Monte Carlo methods [30, l3"Tj . In our simulations we consider two cases - in 

one case A = and the pair potential takes the usual Lebwohl-Lasher (LL) form [33] for 

nematic liquid crystals of perfectly uniaxial molecules which has been extensively studied 

by Zhang et al [34] . In the other case we choose A = 0.2 which represents a biaxial system 

composed of prolate biaxial molecules. For the LL model (uniaxial one) there is a single 

weak first-order N—I transition at a dimensionless temperature T = 1.1232±0.0001 [3~4~[I35]. 

From the Monte Carlo results, as reported in [311 [36], the biaxial model (A = 0.2) exhibits 

a biaxial-uniaxial phase transition at low temperature (T « 0.2) and a uniaxial-isotropic 

transition at higher temperature (T « 1.1). The biaxial nematic-uniaxial nematic transition 

is known to be second order while the uniaxial nematic-isotropic transition is first order. 

The dimensionless temperature has been defined as T = ksTx/e, Tk being the temperature 

measured in Kelvin. 

The interaction of a uniform external magnetic field B chosen along the laboratory Z 
axis (unit vector z), with the ith molecule resulting from its coupling with the longest 
molecular symmetry axis Wi is taken as 

Ur d = -e^{w t .zf-\]. (7) 
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Where, £ is a dimensionless quantity which determines the strength of coupling of the 
molecular symmetry axis with the magnetic field and is given by 

3/j e 

Here, Ak = k\\ — tt± is the anisotropy of the molecular magnetic polarizability and /io is 
the permeability of the free space. In the simulations we take £ to be a positive quantity so 
that the molecules tend to get their long axes aligned along the magnetic field. 
The total energy S of the system is therefore given by 

s = j2 u t p + 2 u l ield ( 8 ) 

(*>i) * 

where, the ( ) bracket represents the nearest neighbours and the two terms on the right 
hand side are obtained from Eqs. Q] and CD 

3 COMPUTATIONAL ASPECTS 

To calculate the thermodynamic observables of interest as a function of A, £ and T we 
have performed a series of Monte Carlo (MC) simulations using the conventional Metropolis 
algorithm on a periodically repeated simple cubic lattice, for four system sizes. A Monte 
Carlo move was attempted by selecting a site at random and then by choosing one of the 
laboratory axes at random and rotating the molecule at that site about the chosen laboratory 
axis following the Barker- Watts method [37]. For generating histograms of energy and the 
constant energy averages of the order parameter and its square, simulations were run in 
cascade, in order of increasing temperature T for a given set of values of A and £. 

We have used 10 6 sweeps or MCS (Monte Carlo steps per site) for the equilibration and 
(4 — 6) x 10 6 MCS for the production run for every set of values of A, £ and T . For the largest 
lattice size (L = 30), the total run length is more than 10 000 times the correlation time. We 
have divided the total run into several (100) blocks by performing independent simulations 
for each set of values of L, A, £ and T so that we could compute the jackknife errors 
[38]. In order to analyze the orientational order we have calculated the second rank order 
parameters (R^n) following the procedure described by Vieillard-Baron [39]. According to 
this, a Q tensor is defined for the molecular axes associated with a reference molecule. For 
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an arbitrary unit vector w, the elements of the Q tensor are defined as 



Qap(w) = ({3w a wp - 5 aP )/2) (9) 

where the average is taken over the configurations and the subscripts a and ft label Cartesian 
components of w with respective to an arbitrary laboratory frame. By diagonalizing the 
matrix one obtains nine eigenvalues and nine eigenvectors which are then recombined to 
give the four order parameters (Rq ), (-R02)) (-^20) an< ^ (-^22) w ith respect to the director 
frame [30] . 

Out of these four second rank order parameters the usual uniaxial order parameter 
(Rqo) (or, S) which measures the alignment of the longest molecular symmetry axis with 
the primary director (n), is involved in our study. The full set of order parameters are 
required to describe completely the biaxial nematic phase of a system of biaxial molecules. 
In our work we have simulated a very short temperature range in the higher sides of T 
within which no biaxial phase occurs. 

We have calculated the reduced specific heat per particle and the ordering susceptibility 
from fluctuations in the energy and the order parameter respectively 

C - {E2) - {E)2 (10) 

x _ N((R 2 00 2 )-(R 2 00 f) (n) 

where E is the scaled total energy of the system. 

Various thermodynamic quantities have been computed using the multiple-histogram 
reweighting method proposed by Ferrenberg and Swendsen [22] . In our simulations we have 
generated the energy histogram for certain values of the dimensionless temperature T and 
corresponding to each energy bin we have generated (constant energy) averages of the order 
parameter and its square. These averages are used to evaluate the order parameter and the 
corresponding susceptibility as a function of temperature using the reweighting method. 
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4 RESULTS AND DISCUSSION 



4.1 Probability distribution functions and free energy for the uni- 
axial and the biaxial models 

In order to analyse the effect of the external field on the phase behaviour of nematic liquid 
crystals we need to simulate the models within a very narrow range of temperature (T = 
1.110 to 1.140) around T/vi- We aimed at getting the transition temperatures with an 
accuracy of ±0.0001 and have therefore used the multiple histogram reweighting technique 
for the purpose. 

We have performed extensive Monte Carlo simulation for different values of the external 
field parameter £ at five or six different temperatures within the said temperature range 
to generate histograms for both the uniaxial and the biaxial models for L = 18, 22, 26 and 
30. For the uniaxial molecules simulations have been performed for five different values of £ 
ranging from to 0.00125 with an increment of 0.000375 for each lattice size, while for the 
biaxial model five different values of £ from to 0.001 with an increment of 0.00025 have 
been used. Thus a total of about 200 simulations were performed for different values of A, 
L, £ and T. The histograms were used to obtain the temperature dependence of internal 
energy (E), order parameter (Rq ) and the corresponding response functions i.e. the specific 
heat (C v ) and the order parameter susceptibility (%). 

The results for the normalized histogram count P(E) = h(E)/ h(E) in absence of 
an external magnetic field for both A = and A = 0.2 have been plotted in Figs. [1] and 
|2] respectively for the L = 30 lattice. These figures show an evidence of double peak like 
structures which are however not well separated and are more so in Fig [2J These merely 
confirm the well known weakly first order nature of the N — I transition which gets weaker 
for the biaxial molecules. With increase in the strength of the magnetic field the overlap in 
the double peak like structures are found to increase even more for both values of A. 

We have derived the relevant part of the free-enegy like functions A(E) from the energy 
distribution functions, P(E) for both A = and A = 0.2 using the relation A(E) = —InP(E). 
For the uniaxial molecules we have shown (Fig. [3]) the field dependence of the free energy 
A at finite size transition for the largest simulated system size (i.e. for L = 30) for which 
the free energy barrier between the two minima is well pronounced and the effect of the 
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external field on AF is clearly visible. Here, the bulk free energy, AF is given by AF(£, L) = 
A(E m ; T, £, L) — A(E±; T, £, L) where E\ is the energy at which either of the two minima of 
A (of equal depth) appears and E m gives the position of the maximum of the free energy 
A. The temperatures in all cases were adjusted to obtain two minima of equal depth. 
We observe that the double-well becomes shallower as £ increases. The external field thus 
reduces the strength of the first order transition and this suggests that for a particular 
value of the field (£ C (L)) the barrier between the two minima of A(E) is likely to vanish 
which would correspond to the end point where the first order phase transition turns into a 
crossover |25j . 

The mesogenic molecules having biaxiality A = 0.2 however does not show any noticeable 
energy barrier separating the two minima in A, even in absence of the external magnetic 
field and for the largest system size (L = 30) simulated (Fig. HJ). This observation thus does 
not provide any conclusive evidence about the nature of the nematic-isotropic transition for 
the system of biaxial molecules. This is true for the system sizes we have investigated and 
it is likely that using systems of significantly larger size one could observe the free energy 
barrier for A = 0.2 and its dependence on the magnetic field. 

4.2 Finite size scaling and the shift in NI transition temperature 
with increasing magnetic field 

For finite systems the peak height of C v increases with increasing system size and the scaling 
relation for C™ ax in a first-order phase transition obeys C™ ax ~ L d [23l [23] . We see that this 
scaling relation holds good for the uniaxial system (Fig. EJ). The data fit well for all values 
of £ used in our simulations. The expected linear fits of ■^ max are presented in Fig. for 
the uniaxial molecules. We have observed that with increasing £ the heights of the maxima 
C™ ax (L) and ^ max (L) are suppressed particularly for the higher lattice sizes L = 26 and 30. 
This observation shows that the presence of magnetic field weakens the first orderedness of 
the N — I transition. 

In case of biaxial system the straight line fits (Figs. [7] and to the data of the peak 
heights of C v and x for £ = 0, 0.00025, 0.0005 and 0.00075 show that the N — I transition 
obeys the first-order scaling laws for this system too. Determining the finite-size transition 
temperature (T^j(L)) from the location of the maximum of the specific heat and the su- 
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ceptibility curves we have used the scaling relation T^j{L) — T^j ~ L~ d and performing a 
linear extrapolation to the thermodynamic limit (L — > oo) the transition temperature T^i 
in the thermodynamic limit was estimated. In Fig. [9] the linear fits of the data for Tni(L) 
obtained from both C v and x plots for the uniaxial molecules are shown and in Fig. [TU] 
the same are plotted for the biaxial molecules. The estimated transition temperatures at 
the thermodynamic limit are listed in Table-1. We find the expected increase in T^i with 
increase in £ for both A = and A = 0.2. The shift STj^i in the transition temperature is 
plotted against £ (£ oc B 2 ) in Fig. [TTJ We obtain good quadratic fits to the data in both 
cases. The quadratic variations of 5T^i with £ do not follow the prediction of the Landau de 
Gennes model for the isotropic-nematic transition which predicts a linear variation of ST^i 
with B 2 . The quadratic functions used for the fits are /(£) = 5.7 x 1(T 6 + 0.57£ + 950.86£ 2 
for A = and g(g) = 8.5 x 1(T 6 + 0.77£ + 1028.57£ 2 for A = 0.2. Therefore, the rise of 
Tni over its zero-field value is steeper for the biaxial molecules than that for the uniaxial 
molecules. For low magnetic fields the slope of the fit for A = 0.2 is 35% higher than that 
for A = 0. The results suggest that the external magnetic field necessary for inducing the 
I(pN) — N phase transition is lower for the biaxial system than the uniaxial system (where 
initial temperatures for both systems are assumed to be equally higher than their respective 
zero-field transition temperatures, T NI (0)). The mean field study of Wojtowicz and Sheng 
[1] for the uniaxial model shows that the slope of a linear fit is lower by 20% in comparison 
to the slope extracted from Fig. 3 of Ref. [7] in which the model under investigation was bi- 
axial. Our observation regarding the shift in T^i with increase in magnetic field is therefore 
consistent with the mean field results. 

The finite size stability limit, T_(L) of the isotropic phase for uniaxial molecules is 
estimated as the temperature where the second, local minimum (at higher energy) of A just 
vanishes as T is gradually lowered below T^i- From extrapolation to the thermodynamic 
limit we estimate T_ for different field strengths as indicated in Fig. [HI by dashed lines. 

4.3 Estimate of the critical field 

An important observation is that the width of the stability limit of the I phase, Tnj(B) — 
T-(B) decreases with increase in magnetic field (Fig. [12]) . From the fit to the data (as 
shown in Fig. IT2"|) we find that the temperature difference T^j(B) — T_(B) varies linearly 

11 



with the field parameter £ (i.e., quadratic in B). With increasing magnetic field the weak 
first order N — I transition of the uniaxial model becomes weaker and T-(B) gradually 
becomes closer to Tnj(B) and at the critical end point they merge into a single point. It 
is, therefore, possible to obtain an estimate of the critical magnetic field by performing 
a linear extrapolation (Fig. [T2]) up to the field strength at which T NI (B) — T_(B) 
becomes zero. For the uniaxial case we get £c = 0.00130. In order to get an estimate of 
£ in real units we may express the magnetic field as B = (3//oe£/ An) 1 / 2 . The energy unit 
e can be estimated by using the experimental and simulated nematic-isotropic transition 
temperatures, i.e. e = )zbTk{B = 0)/Tnj(B = 0). For a common nematic, say MBBA, 
the N — I transition temperature Tk{B = 0) = 320-ft" and the anisotropy of the molecular 
magnetic polarizability A/t w 0.16 x 10 -32 m 3 [UJ. Using these values and the simulated 
value of T NI (B = 0) = 1.1232 for A = we can obtain B « 3040^ T. Therefore, the critical 
magnetic field Be for common nematics corresponding to our estimate, £c = 0.00130, is 
~ 110 T. 

We can also obtain a finite size estimate of £c(L) from the field dependence of AF 
[25] as shown in Fig. [T3j We observe that as the field strength parameter £ is increased 
beyond 0.00125 the free energy barrier decreases much faster. For the uniaxial molecules 
we may estimate the value of £ c by performing a linear extrapolation using the results of 
the barrier height (AF) for three higher values of £ i.e. 0.00125, 0.0013 and 0.0014 (we 
have extended our simulations to £ = 0.0013 and 0.0014 for the largest lattice size L = 30). 
Our estimate of the finite size critical field parameter is ~ 0.00148 which results in a value 
of the critical field Be for the common nematogens to be ~ 120 T. This is larger than 
the estimate in thermodynamic limit given in the previous paragraph. Our estimate of the 
critical magnetic field is of the same order of magnitude (130 T) mentioned in [10] obtained 
from the phenomenological Landau-de Gennes theory and using the experimental values of 
the / parameters for 5CB. But the critical magnetic field predicted by the molecular theory 
of Maier-Saupe [I] is much higher (1000 T) than our estimate. 

The corresponding critical temperature can be estimated from the extrapolation of the 
fitted curve for A = in FigJTTl which gives T c = 1.12555, 5T C being 0.00235. 
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5 CONCLUSIONS 



We have determined the shift in the N — I transition temperature with increasing external 
magnetic field in the thermodynamic limit for two nematic systems composed of uniaxial 
and biaxial molecules. These results represent the first numerical evidence of the first order 
nature of nematic-isotropic transition in an external magnetic field the presence of which 
makes the first order character weaker. The first order nature of the transition is established 
from the finite-size scaling behaviour of the transition temperature T NI {L) obtained from 
the specific heat and susceptibility data. 

The quadratic dependence of 5T NI with £ where £ oc B 2 , for both cases A = and 
A = 0.2, does not agree with the theoretical prediction of LDG model which predicts this to 
be linear for low values of the external magnetic field. The steeper nature of the quadratic 
fit for the biaxial molecules in comparison to that of the uniaxial molecules shows that a 
lower magnetic field strength is capable of inducing the isotropic-nematic transition for the 
biaxial case. 

The free energy like quantity generated by the reweighting technique gives detailed 
insight of the equilibrium properties of the models in presence of an external magnetic 
field. In particular, for the uniaxial case, the effect of an external magnetic field on the free 
energy barrier is clearly observed. Our results show that the width of the stability limit 
of the isotropic phase below T NI (B) decreases quadratically with the external field. We 
estimate the critical magnetic field strength for a system of uniaxial mesogenic molecules 
(~ 110 T). For the biaxial system we have not found any double well structure in the 
free energy. The work further suggests that the first order nature of the nematic-isotropic 
transition of a system of biaxial molecules is considerably weaker than that shown by the 
uniaxial LL model. The presence of an external magnetic field further weakens the first 
orderedness of the NI transition. To study the effect of magnetic field on the nematics 
composed of biaxial molecules systems of much bigger sizes need to be simulated. 
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Figure 1: Normalized histogram for five different temperatures T = 1.121, 1.122, 1.123, 
1.124 and 1.125 for the uniaxial model (A = and L = 30) in absence of an external 
magnetic field. 
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Figure 2: Normalized histogram for five temperatures for the biaxial model (A = 0.2 and 
L = 30) in absence of an external magnetic field. 




Figure 3: Free energy A as a function of energy per molecule with and without external 
magnetic field for the uniaxial model for L = 30. From bottom to top the field parameters 
are 0(B), 0.0009375(D), 0.00125(«) and 0.0013(o). An eighth-order polynomial fit to the 
data is also presented. 
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Figure 4: Free energy A as a function of energy per molecule without any external field for 
the biaxial model for L = 30. 

Table 1: NI transition temperatures (at thermodynamic limit) for different values of the 
external field strength parameter £ for the uniaxial and the biaxial systems. Estimates of 
the super cooling limits T_ are also listed for the uniaxial system. The estimated (jackknife) 
error in each temperature is ±0.0001. 



A 




Tni (from CyvsT) 


Tjvj(from x vs T) 


T_ 







1.1232 


1.1232 


1.1221 




0.0003125 


1.1234 


1.1234 


1.1226 





0.0006250 


1.1238 


1.1239 


1.1233 




0.0009375 


1.1245 


1.1245 


1.1242 




0.0012500 


1.1254 


1.1255 









1.1166 


1.1167 






0.00025 


1.1169 


1.1170 




0.2 


0.00050 


1.1172 


1.1173 






0.00075 


1.1178 


1.1179 






0.00100 


1.1184 


1.1185 
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Figure 5: Finite-size scaling behaviours of the peak height of C v for the uniaxial molecules 
for four different values of the field strength parameter £ = 0.0003125 — 0.00125. The 
estimated errors are of the order of the size of the symbols. 
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Figure 6: Finite-size scaling behaviours of the peak height of x for the uniaxial molecules for 
four different values of the field strength parameter £ = 0.0003125 — 0.00125. The estimated 
errors are of the order of the size of the symbols. 
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Figure 7: Finite-size scaling behaviours of the peak height of C v for the biaxial molecules 
for four different values of the field strength parameter £ = — 0.00075. The estimated 
errors are of the order of the size of the symbols. 
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Figure 8: Finite-size scaling behaviours of the peak height of x for the biaxial molecules for 
four different values of the field strength parameter £ = — 0.00075. The estimated errors 
are of the order of the size of the symbols. 
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Figure 9: Finite-size scaling behaviour shown for two different measures of T NI (L) for the 
uniaxial system at four different values of £ (the data points represented by A and o are 
obtained from the peaks of C v and x respectively). The estimates of T_(L) are also shown 
(□). Extrapolations to the thermodynamic limits are denoted by solid and dashed lines, 
respectively. 
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Figure 10: Finite-size scaling behaviour shown for two different measures of T^j(L) for the 
biaxial system at four different values of £ (the data points represented by A and o are 
obtained from the peaks of C v and x respectively). Extrapolations to the thermodynamic 
limits are denoted by solid lines. 
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Figure 11: Increase in the NI phase transition temperature ST (over zero- field) vs £ for both 
the uniaxial (A = 0) and biaxial (A = 0.2) models. The extrapolation of the fitted curve 
for A = gives ST NI = 0.00235 for the critical parameter £ c = 0.00130 and the symbol (*) 
represents the critical end point. 




Figure 12: Plot of T NI (B) - T-(B) vs £. The solid line is the best linear fit. The estimated 
value of £ c is 0.00130. 
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Figure 13: Free energy barrier height AF(L) vs £ for the uniaxial model and for the lattice 
size L = 30. A linear extrapolation of three nearest points is used to estimate £ c which is 
0.00148. 
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